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' Abstract 

^ (~| Fits to the precision electroweak data that include the NuTeV measurement are 

considered in family universal, anomaly free U{1) extensions of the Standard Model. 
In data sets from which the hadronic asymmetries are excluded, some of the Z' models 
can double the predicted value of the Higgs boson mass, from ~ 60 to ~ 120 GeV, 
^ removing the tension with the LEP II lower bound, while also modestly improving the 

^ confidence level. The effect of the Z' models on both mn and the confidence 

2 level is increased when the NuTeV measurment is included in the fit. Both the original 

0> NuTeV data and a revised estimate by the PDG are considered. 
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To Lev Okun in honor of his 80 'th birthday: a kind and gentle man who does physics as 
he lives his life, with simple honesty and integrity. Although I have told it before, the story of 
my first encounter with Lev bears retelling. It was at the 1976 International Conference on 
High Energy Physics in Tbilisi. Andre Sakharov was not allowed to register for the conference 
but was told that his presence would be tolerated if he wished to attend informally. While 
there was doubtless considerable sympathy for Sakharov among many of the Soviet physicists 
attending the conference, Okun was the only one I saw who dared to associate openly on the 
streets of Tbilisi with Sakharov during the meeting. In another characteristic decision, in the 
period following the fall of the Soviet Union, Lev chose to remain in Moscow to help preserve 
the marvelous school of theoretical physics at ITEP, when he could have accepted offers of 
more comfortable positions outside of Russia. 

Introduction 

When the precision electroweak data began to emerge from LEP, SLC, and Fermilab, 
Okun and collaborators Novikov, Rozanov, and Vysotsky did a complete and independent 
study of the one loop radiative corrections, culminating in the LEPTOP program. [1] Using 
LEPTOP they then made the interesting discovery, contrary to the conventional wisdom of 
the time, that a fourth generation of quarks and leptons is not excluded by the precision 
data, p] They also found that a fourth generation would raise the Higgs mass prediction 
significantly, and that it could resolve the conflict between the SM (Standard Model) predic- 
tion for ntH and the 114 GeV LEP II direct lower limit on mij that arises if the inconsistent 
determinations of the weak interaction mixing angle sin^6'^'^''^ are attributed to underesti- 
mated systematic error in the hadronic asymmetry measurements. |3j In this paper I consider 
a class of Z' models that would also raise the mu predicton into the LEP II allowed region. 
In particular I extend an earlier studyjl] of the EWWG[5] (Electroweak Working Group) set 
of observables, to also include three low energy measurements: the NuTeV vN scattering 
measurement, [B] Moller scattering, and atomic parity violation. [5] Principally as a result 
of the NuTeV measurement, the Z' models have a greater effect on the fits than in the pre- 
vious study: the central value of mn is raised by a factor 2 into the LEP II allowed region 
and the confidence level is modestly improved in some cases. 

The SM Fit 

The SM fit to the precision electroweak data shown in table 1 has two 3cr anomalies 
which together reduce the confidence level of the fit to 2%]^ One is the 3.2cr discrepancy 
between the two most precise measurements of the weak interaction mixing angle sin^^^^^, 
Alr and A\.q, and the corresponding 3.2(T discrepancy between the three hadronic asym- 

^Radiative corrections are from ZFITTER 9 , with two loop corrections to siT?9^^ [TO] and mvi/.[TT] 
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Table 1: SM fits with (A) and without (B) and 
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metry measurements {A^p^, Ap^, Qfb) and the three leptonic ones (A^/j, A^^, 
The second is the measurement of the weak mixing angle, s\.v?6w^'^ scattering by the 
NuTeV experiment. [6 Either or both could be the result of statistical fluctuations, genuine 
new physics, or underestimated systematic uncertainty. If the cause is systematic error in the 
asymmetry measurements, the hadronic asymmetries are by far the more likely candidates 
since they have challenging experimental and QCD-related theoretical issues in common, 
which have no counterparts in the leptonic asymmetry measurements that are measured by 
three quite independent and relatively simple techniques]^ Without prejudging the validity 
of any of the possibilities, in this paper we explore the consequences of the assumption that 
the hadronic asymmetry measurements have underestimated systematic error. We then also 
consider the fit without the hadronic asymmetry measurements, data set B in table 1. 

As observed by Davidson et a/. [13j the NuTeV anomaly is also susceptible to theoretical 
systematic effects from QCDj^In particular, a positive asymmetry in the momentum carried 
by strange versus anti-strange quarks in the nucleon sea would reduce the size of the anomaly. 
Recently the NuTeV collaboration has found a lAa indication for a positive asymmetry. [T^ 
Taken at face value it would reduce the anomaly from 3 to 2(T, as noted by the PDG[15J 
(Particle Data Group), which quotes revised estimates for g\ and g\. In this paper we 
consider both the original NuTeV results and the PDG estimates. 

Table 1, with the original NuTeV data, displays the SM fit with (A) and without (B) 
the hadronic asymmetry measurements, A\.q and A'^pg. In addition to the three low energy 
measurements, we include the usual EWWG data set, with two unimportant exceptions that 
have a negligible effect on the results: the jet charge asymmetry, the least precise of the six 
asymmetry measurements, which we omit because of potential flavor-dependent systematic 
effects that could effect the Z' fits, and the W boson width, which with a 2% error is not 
a precision measurement in the sense of the other measurements that typically have part 
per mil precision. The NuTeV measurement is represented by g\ and which in the SM 
are given by g\ = 1/2 — sin^Ow + 5/9sm'^9w and gj^ = 5/9 sin'^Ow- In fit A of table 1 the 
3.2(7 conflict between Ain and A''pg is manifested by their pulls, +1.8 and —2.7 respectively; 
together with the —2.9 pull of the NuTeV measurement of g1, they are responsible for the 
poor 2% confidence level of the fit. In fit B with A^'p^ and Ap^ excluded, the CL increases 
to 12% but the central value of the Higgs mass decreases to 60 GeV with only 7% CL to be 
in the LEP II allowed region above 114 GeV. 

The best SM fit to data set A using the PDG estimate of g1 and g"^ has /N = 24.2/15 

■^Our fits use the NuTeV measurments of g£ and gj^, which are simple functions of sin^Ow The anomaly 
is manifested in g^, which is more sensitive to sin^Ow- 
•^See [T2] for a more detailed discussion. 

^Davidson et al. also considered an unmixed Z' coupled to i? — 
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implying a still marginal 6% confidence level. The Higgs mass is 85 GeV with an acceptable 
31% CL for rriH > 114 GeV. For data set B the confidence level of the SM fit, shown in table 
5, increases to a robust 35% but the Higgs mass prediction decreases to 58 GeV, with a 6% 
confidence level for mn > 114 GeV. This fit then also provides motivation to consider new 
physics that can increase the mn prediction. 

The Z' Model 

A heavy Z' boson offers a natural solution to an unacceptably light prediction for 
rriH because Z — Z' mixing with a heavy Z' shifts the Z mass downward, corresponding 
to a positive contribution to the p parameter that offsets the negative contribution from 
an increase in the Higgs boson mass. There is a vast literature on Z' models. [16] Here 
we consider the highly constrained class of family-universal models that are anomaly free 
without addition of new fermions to the SM except the e, /i, and r right-handed neutrinos. 
The Abelian generator must then act on SM matter like a linear combination of hypercharge 
and baryon minus lepton number. pTl [T8] which we parameterize by an angle 6x as 



If kinetic mixing occurs [l9l [20] between the Z and Z' bosons we can choose a basis for the 
gauge fields in which it vanishes at the electroweak scale, resulting in a generator of the form 
of equation (1) with a shifted value of 6x — see for instancepT]. Although Y' acts like the 
SM hypercharge Y on SM matter fields, it is a distinct U (1) generator and likewise {B — L)' 
is distinct from the B — L generator of S'O(IO). In particular, we assume the Z' gets a 
large mass from an SM singlet Higgs field with vanishing SM hypercharge but nonvanishing 
Qx charge. This class of models was explored in [18]. Direct experimental bounds were 
extracted from LEP II data in [22] and CDF [23] has obtained even stronger bounds for 
couplings weaker than electroweak, gz' ^ gz/4:- The precision electroweak constraints and 
Higgs boson mass predictions for these models were studied in pi] and ^ for the EWWG 
data set, without the NuTeV, Moller, and APV measurements. 

Gauge invariance of the SM Yukawa interactions requires the SM Higgs boson to have 
vanishing (B — L)' charge, so that Z — Z' mass mixing only occurs through the Y' component 
of Qx- After diagonalizing the mass matrix assuming m|, ^ m| we find[l| the mass 
eigenstates 



Y' 

Qx = cosOx— + sinOx 



{B-LY 
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(1) 



Z 



cos 9m Zq + sin 6m Z'q 



(2) 



cos 6m Z'q - sin 6m Zq. 



(3) 



The mixing angle is 



r cos 6 



'X 



(4) 
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where r and rhz' are the Z' couphng constant and mass scaled to the Z couphng and mass, 
that is 

r = — (5) 
9z 

and 

mz' = . (6) 

mz 

The shift in the Z mass generates a contribution to the obhque parameter T.|25j 

6ml cos^ Ox , . 

otlx = ^ = • (7) 

rriz m^z' 

The effective Zff Lagrangian can then be written as 

^/ = 9z(l + ^) g'fj^f (8) 



where / represents a quark or lepton of chirahty L or R. Since m|/ ^ 1 the mixing angle 
9m is very small, and to leading order in 9m the Zff couphng g'jr is 

g'f=gf + r9Mqic (9) 
where is the Qx charge of fermion /, 

q^x = cos9x y + sin6'x — - — (10) 

and y-^ , , 1^ are respectively the weak hypercharge, baryon number, and lepton number of 
fermion /. The first term in equation (9) has the form of the SM Zff coupling, 

9f = tL - ^^(sin^^^^fi^ + ^oi^sin^^^) (11) 

but with the oblique correction to sin^6'^^^, 

rOB • 2a sm^9wcos'^9w rj. /.o^ 

(5^ sm 9w = ^ r^—aTx. (12) 

cos^o'iy — sm 9w 

For fixed 9x the the effect oi Z — Z' mixing on the electroweak fit is determined by 
a single parameter which we choose to be Tx- The shift in the Zff coupling from the 
Z' admixture, the second term in equation (9), is determined by 

e^r9M = ^. (13) 
cos 9x 
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The fits are obtained by varying Tx in addition to the four SM parameters, mj, ^a^l^{mz)i 
asiiTLz), and mn- The ratio of couphng strength to mass (the effective "Fermi constant") 
is also determined by Tx, 

9z' ^ o'Tx gl , . 

m^, cos^ Ux mz 

or in terms of the scaled coupling and mass, 

aTx 

vcix ^^^^ 



rh\, cos^ Ox 

The LEP II bounds on Z' production [22] can then also be expressed in terms of Tx- Trans- 
lated from the notation of Carena et a/. [22], who extracted the Z' bounds from LEP II 
bounds on contact interactions, the bounds for the first quadrant in Ox are|l] 

- (30.1 + 15^5 tan ^'^^ 
or in terms of the scaled ratio of mass to couphng, 

— > 30.1 cos + 15.5 sin ^x- (17) 
r 

Stronger bounds have been obtained by CDF [23] for Z' couplings that are weaker than 
electroweak as discussed in some specific cases below. 

For the Z-pole measurements the effect of Z' exchange is completely negligible but for 
the three low energy measurements it is of the same order as the oblique corrections and as the 
Z — Z' mixing correction to Z exchange. In addition, the effect on the NuTeV measurements 
includes corrections to the Zvv interaction which is only implicit in the notation. For 
instance, the quoted measurement for the hadronic coupling g\ = gl^+gdi implicitly includes 
a factor 4:gli^ which is equal to one at leading order in the SM. The result for the nonoblique 
corrections, including Z' exchange, is 

^ = ^(1 + tan^x)(sin2^vytan^x - 9 + ISsin^^vK - lOsin^^v^) (18) 
9l 9 

Sgj^ sin^ 



9\ 9 
to which we add the oblique corrections 



[(1 +tan^x)(10 + tan^x - lOsin^^vi/) - 10 - tan ^x], (19) 



= 2aTx(7l - (1 - ^sin^^H^) ■ (20) 
5°^^^ = 2aTx^^ + ^sin^^,^ ■ <5°Wev. (21) 
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There are extensive cancellations among the various Z' corrections to both MoUer scat- 
tering and atomic parity violation, with a simple final result. The sum of the Z' exchange 
correction and the Z — Z' mixing correction to Z exchange cancel exactly with the oblique 
wave function renormalization correction proportional to aTx/'2, from the prefactor in the 
effective Lagrangian equation (8). The only surviving contribution is then the oblique cor- 
rection to sin^6'^''^, equation (12). 

Z' Model Fits 

The Z' fit of data set B with the best confidence level consistent with the LEP 
II limits on rhzi/r is obtained for 6'x = 7r/3. It is compared with the SM fit in table 2 
and figure 1. The decreases by 2.9 units and the confidence level increases from 12% 
to 19%. The improvement can be attributed entirely to the improved fit of the NuTeV g\ 
measurement, for which the contribution decreases by 2.7 units, from 9.3 to 6.6. More 
significant is the improved consistency with the 114 GeV LEP II lower limit on the Higgs 
boson mass. The central value is doubled, from 64 to 126 GeV, completely resolving the 
marginal 7% consistency level of the SM fit with the LEP II limit. The 95% upper limit is 
also doubled, from 124 to 223 GeV. The central value for the ratio of Z' mass to coupling 
constant, rhzi/r = 29, is just consistent with the LEP II lower limit. [22] which from the 
inequality (16) is rhz'/r > 28. For very weak coupling the CDF data|23] for 9x = vr/3 
implies a stronger bound, rfiz'/r > 41, which excludes the fit for gz' < Qz/^ — see table 5 

of m. 

Figure 2 and table 3 show that fits consistent with the LEP II direct bounds on rhz'/r 
are obtained in the entire interval < 6'x < vr/3 while 7r/3<^x<7i"/2is excluded. Table 4 
shows that the fits to the EW data are also acceptable for < < 7r/3, with central values 
for TTiH above 114 GeV and with confidence levels at least as good as the SM. 

The entire Qx parameter space is spanned by the first two quadrants of Ox-, which are 
equivalent to the third and fourth quadrants up to the sign of gz'- In the second quadrant, 
7r/2 < 6*^ < TT, there are no Z' fits which improve on either the SM tjih prediction or the 
confidence level. In particular, for the "x" boson of SO(IO) which couples to Ts/j — (5 — L)/2, 
corresponding to 9x = arctan(— 2), zero mixing is preferred and the x^ is the same as in the 
SM fit. 

Figure 1 and table 4 display "Bayesian" 95% confidence limits for mn which differ from 
the conventionally quoted 95% CL frequentist upper limits that are also shown. Anticipating 
the discovery of a Higgs boson at the LHC, today's Bayesian fit will become the conventional 
frequentist fit when the Higgs boson is discovered. Since we will then not scan on uih there 
is one additional degree of freedom. The Bayesian 95% interval is the domain of mn with 
confidence level > 5% for fits with inn fixed and indicates the allowed range in mu after the 
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Figure 1: distributions as a function of mn for the SM fit (dashed hues) and the Z' (sohd 
hues) fit with Ox — 7r/3 as shown in table 2. For each fit the lower horizontal line is 
the symmetric 90% confidence interval and the upper horizontal line is the Bayesian 95% 
confidence interval defined in the text. The vertical dash-dot line indicates the LEP II direct 
lower limit on mu- 
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Table 2: The SM fit of data set (B) compared to the Z' fit with = tt 
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Figure 2: irhzi/r as a function of 9x for the Z' fits (solid line) compared with the LEP II 
95% lower limits given by equation (17). 



Model 


rhzi /r 


rhz'/r [LEP II] 


r=i(TeV) 


^x = 


50 


> 30 


4.6 


dx = 7r/12 


48 


> 33 


4.4 


Ox = 7r/6 


45 


> 34 


4.1 


Ox = 7r/4 


37 


> 32 


3.4 


Ox = vr/3 


29 


> 28 


2.6 



Table 3: The value of rhz'/r from the electroweak fit compared with the LEP II lower bound. 
The last column shows the Z' mass if gz' = gz, i-e., if r = 1. 
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Higgs boson is discovered. In contrast, the usual frequentist 95% upper limit (the maximum 
of the symmetric 90% confidence interval) uses the best fit with mn scanned as a free 
parameter to estmate the likelihood that the Higgs boson will he discovered within the given 
range of values. Table 4 shows that the Bayesian limits lie above the frequentist limits and 
that they are also increased in the Z' models relative to the SM. 

The Z' fits to data set A, which includes the hadronic asymmetries A\.q and A^^, are 
not shown, because they are little changed with respect to the SM. For most of the Ox 
domain they are not changed at all. In the most favorable fit only decreases from 28.4 to 
28.0, with a confidence level of 0.014, less than the 0.019 confidence level of the SM fit, and 
with only a modest increase in the Higgs boson mass, from 94 to 130 GeV. 

As observed in |13j the results presented by NuTeV [6] for sin^^^y and g\ ^ are sensitive to 
the assumption that the nucleon strange quark sea is symmetric between s and s quarks. The 
NuTeV collaboration has since found a 1.4cr indication of an asymmetry in the momentum 
carried by strange quarks, 5*" = +0.0020(14). [H] The collaboration has not revised their 
value for sin^Ow based on this result, but they do state that an asymmetry of +0.007 would 
make their measurment of siia^Ow consistent with SM prediction. Using this statement and 
taking the new S~ measurement at face value, the PDG[T2] has estimated revised values for 
siri^ 9w, gl, and gj^, which we use in the alternative fits shown in figure 3 and tables 5 and 
6|^ With these estimates the disagreement with the SM is reduced to a 2a effect. 

As can be seen by comparing the SM fits of tables 1 and 5, the effect of the revised 
estimate of the NuTeV results on the global fit to data set B is to raise the confidence level 
from 0.12 to 0.35, while the inconsistency with the LEP II constraint on ttih remains: the 
central value is tuh = 58 GeV and the confidence level for tuh above 114 GeV is 6%. Again 
the Z' models can resolve the inconsistency with the LEP II lower bound on tuh while 
maintaining the confidence level of the fit. In particular, for 6x = vr/3, shown in figure 3 
and table 5, the central value doubles, to 109 GeV, and the CL for niH > 114 GeV increases 
to 45%. All Z' models in the domain < < ''r/3 increase the rriH prediction into the 
LEP II allowed region, with acceptable fits to the precision data and with values of mz' / gz' 
consistent with the LEP II lower bounds as shown in tables 6 and 7. 

For data set A with the PDG estimates of gl and gj^, the SM fit yields xV^ = 24.2/15 
with a still marginal confidence level of 6%. The Higgs mass prediction is acceptable, at 
ruH = 85 GeV, with 30% likelihood for ttih > 114 GeV. In most cases the Z' fits prefer zero 
mixing, and in all cases the effect on and mn is negligible. 



^The anomaly in sin^ew is 0.2277 - 0.2227= 0.0050 so the shift from S' = +0.0020(14) is 2/7 of the 
anomaly or —0.0014(10). Combining errors in quadrature then gives the values quoted by the PDG. 
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Tabic 4: Best fit value of Tx , confidence level, and mH predictions for SM and Z' fits to 
data set B. 
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Figure 3: distributions as a function of itlh for the SM (dashed lines) and Z' (solid lines) 
fits shown in table 5, with PDG values for g\ and g\. For each fit the lower horizontal line 
is the symmetric 90% confidence interval and the upper horizontal line is the Baycsian 95% 
confidence interval defined in the text. The vertical dash-dot line indicates the LEP II direct 
lower limit on m^/. 
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1.3 


Rb 


0.21629 (66) 


0.21584 


0.7 


0.21581 


0.7 


Rc 


0.1721 (30) 


0.1722 


-0.05 


0.1723 


-0.07 




0.923 (20) 


0.935 


-0.6 


0.935 


-0.6 


Ac 


0.670 (27) 


0.669 


0.04 


0.669 


0.03 


9l 


0.3010 (15) 


0.3043 


-2.2 


0.3037 


1.8 


9l 


0.03080 (11) 


0.03003 


0.7 


0.03020 


0.5 




0.23339 (140) 


0.23117 


1.6 


0.23144 


1.4 


xwiCs) 


0.22939 (190) 


0.23117 


-0.9 


0.23144 


-1.1 


rrit 


172.6 (1.4) 


172.3 


0.2 


172.3 


0.2 




0.02758 (35) 


0.02754 


0.1 


0.2768 


-0.3 


as{mz) 


0.118 




0.120 






rriH 




58 




109 




CL{mH > 114) 




0.06 




0.45 








118 




201 




Tx 








0.033 




rhz'/r 








31 




xVdof 




14.4/13 




12.4/12 








0.35 




0.42 





5: SM and Z' {Ox = 7r/3) fits to data set B with PDG values for gl and 
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Model 


Tx 




CL(x^) 


rriH 


CL{mH > 114) 


m^'^"(Frcq.) 


m^'^"(Bayes) 


SM 




14.3 


0.35 


58 


0.06 


118 


192 


^x = 


0.043 


13.7 


0.32 


104 


0.42 


189 


285 


9x = vr/12 


0.044 


13.3 


0.35 


109 


0.42 


197 


300 


9x = n/6 


0.043 


13.0 


0.37 


114 


0.50 


203 


312 


9x = 7r/4 


0.039 


12.7 


0.39 


114 


0.50 


202 


317 


9x = 7r/3 


0.033 


12.4 


0.42 


109 


0.45 


201 


320 



Table 6: Best fit value of Tx , confidence level, and mn predictions from SM and Z' fits 
to data set B with the PDG estimates for gj^ and gj^. 



Model 


rhz'/r 


mz'/r [LEP II] 


"■iz'|r=i(TeV) 


9x = 


55 


> 30 


5.0 


9x = 7r/12 


52 


> 33 


4.7 


9x = 7r/6 


47 


> 34 


4.3 


9x = 7r/4 


40 


> 32 


3.6 


9x = 7r/3 


31 


> 28 


2.8 



Table 7: The value of rhz'/r from the fit with PDG estimates for g\j^ compared with the 
LEP II lower bound. The last column shows the Z' mass if gzi = gz, i-e., if r = 1. 
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Conclusion 

We have extended a previous study of a simple class of Z' models to include three low 
energy measurements in addition to the EWWG set of Z-pole observables. In the earlier 
study we found that the Z' models had little effect on the o^' oii the central value of 
rriH but that they did significantly extend the upper limit on rriH for the data set without 
and Ap^, increasing CL{mH > 114GeV) from 2% in the SM to as much as 30% in 
the Z' models. With the inclusion of the low energy measurements and, in particular, the 
NuTeV data, the effect of the Z' models is more pronounced. Both the central value and 
the upper limit of the Higgs boson mass are increased significantly, by as much as a factor 
two, and the confidence levels are modestly improved. If underestimated systematic error 
proves to be the explanation of the 3.2cr difference in the SM determination of sin^^^'^'^ from 
the leptonic and hadronic asymmetry measurements, then the Z' models discussed here can 
resolve the resulting problematic prediction for the Higgs boson mass, providing completely 
satisfactory fits to the precision EW data. Though it remains to be confirmed by a careful 
study, we expect that the Z' models favored by the fits can be excluded or confirmed at the 
LHC. 

This work was supported in part by the Director, Office of Science, Office of High Energy and 
Nuclear Physics, Division of High Energy Physics, of the U.S. Department of Energy under Contract 
DE-AC02-05CH11231 
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